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["t I The Rattleback is a very popular science toy shown to students all over the world to 

demonstrate the non-triviality of rotational motion. When spun on a horizontal table, 
this boat-shaped object behaves in a peculiar way. Although the object appears sym- 
metric, the dynamics of its motion seem very asymmetric. When spun in the preferred 
4^ direction, it spins smoothly, whereas in the other direction it starts to oscillate wildly. 

^ The oscillation soon dies out and the rattleback starts to spin in the preferred way. We 

^ will present an analytical model capable of explaining this behaviour in a simple and 

intelligible way. Although we aim at a comprehensive treatise, we will study the details 

Q only when they are necessary to understand the calculation. After presenting the calcu- 

lations we will discuss the physical validity of our assumptions and take a look at more 
sophisticated models requiring numerical analysis. We will then improve our model by 

^ assuming a simple friction force. 

Qh Keywords: Rattleback, celt, Celtic stone, unidirectional, asymmetry 

> 1 Background 

For centuries the mysterious properties of the rattleback were used to predict the future 
\^ and bring messages from other worlds. The alternative name "Celt" comes from druids, 

^ or Celtic priests, who are thought to have used the rattleback in their rites. This was 

probably a very good way to convince your followers of your powers, as the surprising 
behaviour looks quite magical even in our times. The priest was also able to adjust 
. • the predictions according to his needs by spinning the rock in the right direction. The 

peculiar motion of the object is still sometimes associated with all kinds of strange fields 
and ether forces. More conventional wrong explanations including the Coriolis force and 
^ the tennis racket theorem are also common, but fail to describe the situation correctly. 

*This article is based on a bachelor's thesis with the same title 



1 



The first mathematically satisfactory analysis was written by Sir Gilbert Walker. In 
his paper "On a dynamical top" (1896) he assumed the celt to roll without slipping and 
neglected all dissipative forces. Walker also examined the motion assuming purely oscil- 
latory initial values. Walker's paper was not available to the author, but other sources 
(such as Garcia&Hubbard) mention the analysis to be quite similar to the modern trea- 
tises with the same assumptions. Walker's paper is criticized for being quite difficult to 
follow, but his analysis has proved to be very resilient. After the seminal paper of Walker, 
some celt-related papers were published in mostly educational publications. Interest in 
celt dynamics began to grow in the seventies and the 1980's became the golden age of 
celt dynamics. Many of the most important papers on the rattleback were published 
during a relatively short time. Garcia&Hubbard trace this sudden increase in interest 
to the appearance of J. Walker's popular article in Scientific American. These articles 
concentrated mainly on numerical simulations of different celt models with dissipative 
forces and their significance in the dynamics. 

In 1986 Sir Hermann Bondi published his classic paper "The Rigid Body Dynamics 
of Unidirectional Spin", in which he proposed an improved analytical model for the 
rattleback. Bondi's model forms the basis of most analytical treatises, although some of 
the credit should probably be given to G.T. Walker. Bondi mentions the paper of Walker 
to be the only treatise he has been able to find, despite several papers were published in 
the years between these two analyses. Bondi does not properly investigate the validity 
of his assumptions and the motivation of the spin reversal is not very rigorous. These 
are better addressed by Garcia&Hubbard, who study the validity of the assumptions 
numerically and propose a better model for the spin reversal. Some other papers have 
also used numerics to investigate the physical requirements of the conservative model. 
We will take a look at these investigations after our main treatise. 

Some alternative proposals are quite unrealistic as such: Caughey's celt model men- 
tioned by Blackowiak&al and Garcia&Hubbard combined a viscous drag with springs 
parallel to the coordinate axes and, according to Caughey, could "capture the essential 
features of the rattleback" . However, Caughey's model could describe the motion per- 
fectly and still fail to explain the behaviour of the celt, as the system in question simply 
is not a celt with any reasonable assumptions. 

Modern research in celt dynamics has become quite theoretical, as can be seen from 
the paper by DuUin&Tsygvintsev whose treatise could well be regarded as mathematical 
physics. Borisov&Mamaev have also used modern techniques in their paper "Strange 
attractors in rattleback dynamics". With certain energies, the celt is shown to be an 
example of classical chaos in a conservative system. Despite its title, the paper also 
investigates the overall behaviour of the celt by using symmetries and phase mappings. 
The variables used in this paper are not very easy to link with the physical situation, 
which makes the interpretation of the various diagrams quite difficult. 

One of the biggest strengths of the model presented here is the choice of simple vari- 
ables. We will follow the contact point between the celt and the table, which can be 
expressed in the coordinates fixed to the celt. This choice is much more convenient than 
the usual Euler angles(e.g. Markeev), not to mention the Andoyaer-Deprit -coordinates 
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of Borisov and Mamaev. The fact that these kind of papers are stiU being published 
clearly shows that the classical mechanics of realistic systems is far from being completely 
sorted out. 



2 Equations of motion and unidirectionality 

We will now present a simple model for the celt behaviour. We shall begin by looking at 
the stability of rotation and continue with an improved analysis of the spin reversal. To 
explain the observed unidirectionality, we will present an adaptation of Bondi's classic 
treatise which forms the basis of most studies of celt motion. We shall try to add some 
comments to the equation-dominated presentation of the original article. 
Let us first introduce the notation: 



s = position of the center of mass 

r = vector from the center of mass to the point of contact 
F = force exerted to the celt by the table 
Lo = angular velocity of the celt 
h = angular momentum of the celt 
M = mass of the celt 

u = upward pointing unit vector perpendicular to the table 
Using these we can write the equations of motion 

= F - Mgu (1) 
v + a;xr = (2) 

which give 

-=Mrx[-+,u]. (4) 
Taking a scalar product with u and swapping the operations then leads to 

d , s , . dv , ^ 

-(u.h) = M(uxr).-. (5) 

At this point Bondi discusses the interpretation of this equation in ordinary rotational 
motion. We shall move straight into the rattleback case and define the form of the 
rattleback accordingly. Wc shall thus define an object whose axes of inertia differ from 
the symmetry axes. Bondi arranges this by fixing the coordinates to the center of mass 
and choosing the axes of inertia to be the coordinate axes with the z-axis pointing 
downward at rest position. The moments of inertia are A, B, and C with A > B. The 
form of the object is then defined to be 



i--p(-Y-Q"4-l;s(yy 



(6) 
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where a is the equiUbrium distance form the center of mass to the point of contact. This 
form can also be taken to represent the form only to second order near the equilibrium, 
in which case the object does not have to be an ellipsoid. First order terms are absent, 
as the point of contact is an extremum. 

Physical validity limits the choice of the parameter values. The surface has to be 
concave, which requires 

p > 0, s > and ps > q^. (7) 

This corresponds to a maximum of z at the origin. 

We must also require the equilibrium position to be stable for the oscillation to take 
place. This means that the radii of curvature have to be larger than the equilibrium 
distance a, which is satisfied if 

1 > p, 1 > s and (1 -p){l - s) > . (8) 

If the second order form is only an approximation, there has to be a safety margin 
to ensure stability during the peak amplitude. For typical celts both requirements are 
obviously satisfied. We shall thus assume the parameters to be physical and reasonable. 

The celt will be analyzed by using a system of coordinates fixed to the celt itself. The 
contact point is now given by the vector r = (x, y, z) where z is obtained from ([g]). The 
unit normal vector of the table can be obtained from the definition, as the unit normal 
of the celt surface at the point of contact has to coincide with the unit normal of the 
table. Taking the gradient of ^ results in 

(px + qy qx + sy \ 
u = -w[ , ,1 , (9) 



where w the normalization factor given by 



= y ('^)' + ('^)' + 1- (10) 

The normal vector and the third coordinate axis point at opposite directions at equilib- 
rium, as expected. 
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As u is constant in inertial coordinates, we can write 

= ^ = u + uxu. (11) 

Taking a cross product with u and using the expansion formula for triple product results 
in the expression 

(J = ii X u + nu, (12) 

where n is the "spin" or the vertical component of the angular momentum. We can 
expand the components of the angular momentum assuming x and y to be small enough 
for higher order terms to be dropped. The spin is assumed to be reasonable but it is not 
explicitly assumed to be small. 
The z-component dominates if it is present, which enables us to write 

wi = - [qx + sy - n(px + qy)] 
a 

W2 = - [-px - qy - niqx + sy)] (13) 
a 



-^{ps - q ){xy - yx) - n[l - -{ ) -- ) 

2 a 2 a 



The coordinate axes are the axes of inertia, which allows us to write 

hi = Aoji, /i2 = Buj2, and = Cu^. (14) 
From Q we obtain at first order 

vi = [px + qy + n{qx - (1 - s)y)] 

V2 = [qx + sy + n{{l - p)x + qy)] (15) 
^3 = second order. 

Here all the terms arise from the zeroth order terms in z and the third component of the 
angular momentum. By inserting these expressions to Q and taking into account that 
dP/dt = P + a; X P for all vectors P, we can write the first component of the equation 

as 

Aui — {B — C)(jJ2U)3 
A Tl 

= — [qx + sy- nipx + qy) - hipx + qy)] iB-C)\px + qy + niqx + sy)] . . 

a a (Id) 

= —Ma{qx + sij + n(l — 2p)x — 2nqy + n[(l — p)x — qy] — v?[qx — (1 — s)y\} 

+Mga[qx — (1 — s)y] 

Here the first equality comes from expanding the expression to first order. The first term 
is the derivative of angular momentum and the second term comes again purely from 
the the third component of angular velocity. Second equality is the equality of equation 
Q. When expanding the expressions, it is convenient to remember at each step that we 
are only keeping terms up to first order. 
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Similarly we obtain the second and third components of Q 

BUJ2 — (C — A)uJ3ljJl 

B n 

= [px + qy + niqx + sij) + hiqx + sy)] + - iC - A)\qx + sij - nipx + qy)] 

a a [u) 

= Ma{px + qij + n^qx — (1 — s)y\ + h[qx — (1 — s)y] + n2[(l — p)x — qy]} 

-\-Mga[x{l — p) — qy] 

and 

Cuj3 -{A- B)ujiuj2 = -Cn[l - -( - - ) 

2 a la 

1 1 

H — T^CnUpx + qy)(px + qy) + (qx + sy)(qx + sy)] + -^Cips - q^){xy - xy) 

= Mx{qx + sij + n[(l — p)x — qy]} — My{px + qij + n{qx + sy) + h[qx — (1 — s)y]} 
+Mn{x[{l — 2p)x — 2qy] — y[2qx — (1 — 2s)?/] — xn[qx — (1 — s)y] — yn[{l — p)x — qy]} 

-m-[q{x'^ + y^) + (p - s)xy]. 

(18) 

The coefficient of h is of zeroth order, whereas the rest of the equation is of second order. 
This means that h itself must be of second order, which can be used to simplify the 
other two equations. We will further simplify them by adopting the following shorthand 
notation 

A + Ma^ = aMa'^,B + Ma^ = pMa^, C = -fMa^. (19) 



These allow us to write the equations ( |16[ ) and (17) in the form 

a{qx + sy) — {a + (3 — j)n{px + qy) — {(3 — 'y)n'^{qx + sy) + nx + n'^y = - [qx — {l — s)y] (20) 

a 

and 

P{px + qij) + {a + P — j)n{qx + sy) — {a — 'j)n'^{px + qy)—ny + n'^x = -[qy—{l—p)x] (21) 

a 



We will need equation ( 18 ) only with small values of n, when it can be written as 

a^jh = j(ps - q^){xy - ijx) + [qiyij - xx) +pyx - sijx] + -[(p- s)xy + g(x^ + y^)] (22) 

a 



Now we can take a look at the solutions of the simplified equations of motion (20) 
and(21). As we are examining an oscillation around the equilibrium point, it is natural 
to assume a solution of the form Ce"^. By inserting this solution to the equations and 
cancelling the exponent, we arrive at a system whose characteristic equation is 

((T^ + n^)[a/3(ps - q^)a'^ + (a - fi)qncr] 

+ (cj2 + n^){n^[l - p{a - 7) - - 7) + {ps - q^){a - l){P - 7)] + ^[as + /3p - {ps - q^){a + 13)]} 
2 

+ ^[2 - (1 + a + /3 - 7)(P + s) + 2(a + /? - 7)(ps - q^)] + {^-)\l - p){l - s) - q^] = 
a a 

(23) 
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This equation determines whether the system has nonzero solutions. The character- 
istic is of fourth order, which would allow for a explicit solution by using the Cardano 
formulae. Fortunately, we do not have to solve this equation to understand the behaviour 

of the celt as the important features can be seen from the equation itself. The roots of 
the equation depend on the direction of spin in a straightforward way, namely: 

If a is the solution of the characteristic equation for some value of n, 
the solution for the value -n is -a. 

This is easy to see, as the variables are always squared except in one term in which 
they are both present. It is worth noticing that our model is fully conservative. The 
diminishing amplitude is not caused by dissipation, but is an indication of energy flow 
between the oscillation and the rotation. 

The behaviour depends on the overall distribution of signs. The main classes are the 
following: 

(a) The real parts of the roots all have the same same sign on some interval of spin 
values. This makes one of the spin directions stable and the other unstable. 

(b) The real parts of the roots have different signs on some interval of spin values, 
which leads into unstable motion in both directions. In both directions energy flows 
into the oscillations corresponding to positive real parts. Some positive roots may have 
much smaller absolute values than their negative counterparts, which would result in 
the rotation momentarily gaining energy from the oscillation. 

We are mostly interested in case (a), which essentially the explains the behaviour of 
commercial rattlebacks. The direction of energy flow depends on the direction of spin, 
which causes the observed difference in stability. 

According to Bondi's linearized model, unidirectionality is observed when a body of 
suitable form has a stable equilibrium and unequal principal moments of inertia with axes 
of inertia differing from axes of symmetry. If these conditions are not met, the coefficient 
(a — multiplying the important cross term will vanish. In fact, unidirectionality 
requires the radii of curvature to be different, too. All the requirements are beautifully 
expressed in Bondi's further studies and in the analysis of torque by Garcia&Hubbard. 
We will study both of these later on in this paper. 

The difference in the axis alignment is caused by the coefficient q, but the object can 
also be unsymmetrized by adding weights to a symmetric body. Most experiments have 
used commercial rattlebacks, whose properties have to be estimated from the dynamics 
itself. Taking a symmetric body and adding weights into it would provide us with an 
object with known inertia tensor and shape. This would enable us to make more accurate 
comparisons between the simulations and observed behaviour. The author has made a 
few celts this way and they work surprisingly well. The post-reversal spin rates are often 
left smaller than in commercial celts but the unidirectional behaviour is quite easy to 
observe even with relatively rustic models. 

Bondi then analyzes the relation between the sign distribution and celt shape. This 
results in a diagram classifying celts into types 0,1 and II, although the zero region 
was not given a special name in the original article. In this classification 0-type cor- 
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responds to instability in both directions and I-type to unidirectionality. According to 
Garcia&Hubbard, type II behaves similarly to the unidirectional type, but the axis of 
rotation differs from the vertical equilibrium position. The behaviour of the celt thus 
depends on its position in the chart. The parameters of the celt affect its stability and 
by fixing some of them we can find threshold values for the rest. In addition to the 
linearized model, this belongs to the classic material of Bondi's paper. Some aspects of 
this analysis are presented in the appendix. 

After classifying the various celt types, Bondi considers the spin reversal. Rotation 
in the unstable direction comes to a halt while an oscillation remains. Equation (22) 
shows, that the spin cannot remain zero, but the celt starts to spin again. According to 
Bondi, it is very implausible that the object would start spinning in the same direction 
after having just stopped due to the instability. The object would thus start spinning 
in the stable direction and the negative real parts would pump the energy into the spin. 
In principle, Bondi's model thus explains the observed spin bias. 



3 Torque and the spin reversal 

The explanation given in Bondi's final conclusions is quite credible but as an essential 
feature of celt dynamics the spin reversal deserves to be studied more thoroughly. The 
reversal is investigated in detail by the mostly numerical paper by Garcia&Hubbard, in 
which the authors present in a rather concise way a model based on average torque. We 
will now go through their reasoning and comment on its details after the calculation. 

Similarly to the previous model, Garcia&Hubbard neglect dissipation and assume the 
celt to be rolling without slipping. The initial situation is a slowly rotating celt with 
energy stored in the oscillation. 

We are interested in the vertical torque which is obviously equal to the vertical com- 
ponent of the derivative of angular momentum: 

= u • dh/dt = Mu • (r x dw/dt) = M(u x r) • dw/dt (24) 

Integration by parts gives 

j u - dh = = [(u X r) • v]r - y V • (u X dr + du X r) (25) 

As the motion is approximately periodic during the interval T, we obtain 

[(u X r) • v]r = (26) 
By writing v = r x and r = Ru we can express this as 
T 



/ ((Ru) X w) • (u X (Rdu) + du X (Ru)), (27) 
Jt 
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where the matrix R has entries 

1 



rii/a 



[(1 + ^l:)/9 + (1 - ^)/</.] = - cos' i/e - sin' 



ri2/a = rsi/a = -^(l/<^ - 1/^)(1 - V')^ = "(V^ " V^) cos^sine 
r22/a = [(1 - ^)/0 + (1 + ^)/</'] = - sin2 i/O - cos^ i/(t> 



(28) 



J^lS = -^23 = ^^31 = -^32 = 
^^33/0 = -1 

The original paper lacks a minus sign in the first equation. This is probably just a 
printing error, as all the results are consistent with the correct form. 

Garcia and Hubbard define the curvature parameters as follows: 
9 = the ratio of a to the smaller radius of curvature 
(j) = the ratio of a to the larger radius of curvature 

= cos^ ^ — sin^ ^, 

where is the small angle between the direction of the smaller radius of curvature and 
the axis of inertia deviating from it. 

Bondi uses the the same parameters but defines them in a completely different way. 
The definitions are equivalent, as can be seen by diagonalizing the quadratic form and 
using the formula for curvature radius. The cosine of the skew angle can be computed 
from these results and the use of elementary trigonometry gives the result. Some of 
Bondi's results have to be reformulated to see the connection better. We did not present 
the classification of celt types explicitly but the equivalence of the definitions and the 
calculations involved are quite illuminating as such. Among other things, they explain 
the stability condition Q postulated by Bondi. Despite its algebraic nature, we thus 
present the analysis in the appendix. 

Expanding the integral up to fourth order results in the expression 



T 

M 



Tv = a [uJidui{a + rn) + ujidu2ri2 + uJ2duir2i + U2du2{a + r22)] 
JT 



(29) 



where the quantities involved are obtained from linearized equations. 

The derivative of the state vector [wi uj2 ui U2] can be obtained from the vector itself 
by multiplying with the matrix 



(30) 









d 


e 








f 


k 





1 








-1 












where the coefficients are defined as 

g{a + r22) 



d 



gru 



g{a + m 



gru 
Pa? 



(31) 
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The eigenvalues of A are 



^1,2 



-{(e-/)±[(e-/)2 + 4(e/-dA;)]2} 



(32) 



If the geometric axes coincide with the axes of inertia, we obtain the squared frequencies 
e and — /. We thus define the two eigenvalues to be 2 = 1 /'• The initial conditions 
are obtained from the eigenvectors: 



V 



-e{e' + f) + dk -e{e' + f) + dk f'd f'd 

-e'k -e'k f{f' + e) - dk f{f + e) - dk 

i^fe>k -iVe'k ij^if + e) + e) 



Initial conditions 2(^1 + V2) and2(v3 + V4) give the torques 

= jMa2a(i - - i)(l - i:^)he'Me' + f) - dk]vl 
f,f = \Ma'p{^ - - i)(l - + e) - dk]v}, 



(33) 



(34) 



where vi_2,3,4 are the column vectors of V and Vf,Ve are normalization factors with 
dimension s^. 



4 Commentary 

The first part of the calculation by Garcia&Hubbard is more or less a straightforward 



application of rigid body dynamics. Expression (24) has the same form in the rotating 
frame of reference, which can be shown by using simple vector algebra. Matrix R can 
be derived from Bondi's equation ^ by inverting the relation and using equations (36) 



to transform the result to these variables. The calculation becomes easier, if one notices 
the relation ps — = 6(j) . The entries of the dynamical matrix A can be derived from 
([3]) by using the relation between r and u. The horizontal forces near the equilibrium 
point can be found by using the components of u, which results in the expression given 
in Garcia&Hubbard. 

After obtaining the various components, the calculation becomes rather mechanical. If 
the skew angle vanishes, the entries of the matrix can be deduced by drawing a picture. 
The torques can be found quite easily and the result can be written down. In both 
methods one has to pay special attention to the signs, especially as the vertical axis 
points downwards. The skew angle enters the expressions in A and R only because the 
coordinate axes differ from the curvature axes. In most celts this angle is quite small, 
which allows us to study the connections between the various quantities more easily by 
temporarily assuming = 0. This assumption is not possible in our actual treatise, as 
the asymmetry is necessary for the celt to work. 
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The substitution term in the integral is assumed to vanish due to the periodicity. The 
motion is not absolutely periodic and the forces considered rather small. The motion is 
approximately periodic during one cycle, which makes the substitution negligible. 

The linearized equations mentioned by Garcia&Hubbard are ambiguous, as their paper 
does not explicitly present clear candidates. One possibility is to use equations ([9]), (10), 
and ( 13 ) assuming the spin to be small. Taking terms up to second order into account one 
obtains the same form of equation (29). It can be shown that the squared frequencies -f 
and e' are positive real numbers. The choice of signs ensures that the primed frequencies 
take their unprimed values when the skew angle is taken to zero. The eigenvectors are 
obtained from A and choosing a suitable linear combination results in a trigonometric 
form in both modes. According to the dynamical matrix the derivatives of u are the 
angular frequencies, which can be verified from the eigenvectors. This connection comes 
handy in simplifying the expressions. The normalization factors fix the amplitude of the 
oscillations and correct the dimension of the initial values. The diminishing amplitude 
of the oscillation is taken into account through the dynamics included in the linearized 
equations, which causes the final state to differ from the initial one. The averaging leads 
into integrating the squares of trigonometric functions over the period and a nonzero 
result. Simplifying the result after the integration is quite tedious, but results in the 
expression obtained by Garcia&Hubbard. We have presented this analysis mainly to 
verify Bondi's conclusions, which makes the technical details and the exact form of the 
result not so important. The calculation verifies that a suitably asymmetric body with 
two unequal principal moments of inertia and two differing radii of curvature oscillating 
on a horizontal table experiences a nonzero torque, which causes the spin reversal. It 
is very pleasing to note that exactly the same conditions are obtained by Bondi during 
the classification of celt types. The intermediate steps presented in the commentary 
probably differ somewhat from the original path taken by Garcia&Hubbard. We could 
thus claim the conditions to be obtainable in three at least somewhat different ways. 
The diminishing spin and increasing oscillation are caused by the torque, which remains 
nonzero even after the celt comes to a halt. This nonzero torque is responsible the spin 
reversal. In the stable direction energy flows back to the rotation, which can already be 
seen from the preceding model by Bondi. 



5 Validity of models 

With the assumptions made here, Bondi's linearized model by seems to be a generally 
accepted explanation for the unidirectional stability of the celt. On the other hand, the 
true nature of the spin reversal has been discussed more frequently and many numerical 
papers give a larger role to the dissipative forces ignored by Bondi. Garcia&Hubbard 
use numerical simulations to show the negligibility of aerodynamic drag in the motion 
of typical celts made of relatively dens e materials. Although Garcia&Hubbard use a 
relatively coarse linearized model, we can rule out the effects of the airfoil form. Another 
class of models allow the contact point to slide. Studying these models is probably 
in order, since maintaining a non-sliding contact seems to require unrealistically high 
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coefficients of friction. The effects of sliding friction are investigated in the papers by 
Magnus, which seem to be unavailable. Taking into account the sliding friction and 
sliding in all directions will naturally lead into considerable algebra. It is therefore 
not very surprising that Lindberg&Longman report the discussion in Magnus to be 
too heavy to offer any insight into the motion. Lindberg&Longman's simulations with 
the assumptions made in Bondi showed unidirectionality, which confirms the physical 
validity of the model. 

The coefficient of friction can be taken to be very large, if wc model the point of 
contact to be a mathematical point. This would prevent the slipping without slowing 
the celt down. This is not necessary or even desirable, as realistic celts do slow down. 
The various celt models also seem to predict extra reversals, if the simulation is run 
for unrealistically long times. In reality, celts typically reverse only once before coming 
to a halt due to friction. In the paper by Lindberg&Longman, the extra reversals take 
place at 125 seconds and 150 seconds, depending on the initial direction of spin. If the 
graphs are restricted to realistic spin times,however, we see a beautiful reversal followed 
by rotation in the stable direction. Adding a dissipative vertical torque without any 
other dynamical effects would thus result in a rather realistic model. This could be done 
by replacing the mathematical contact point with an area of contact, which is in fact 
a better description of a contact between two objects. The celts built by the author 
seem to work the best when the initial spin is rather small and the oscillatory motion 
demanding the largest coefficients of friction thus not so violent. On the other hand, 
the existence of a vertical torque caused by the friction is clearly visible. The modest 
oscillation left behind cannot overcome the friction and the celt is left non-rotating. 
Although multiple reversals do not occur in our linearized model, assuming an area of 
contact still makes the model more realistic. 

We can conclude that our assumptions are reasonably physical and our treatise thus 
a satisfying explanation for the dynamics of the celt. According to simulations, higher 
order calculations would result in quantitatively reasonable results. Taking higher orders 
into account would make the necessary analytical calculations rather tedious. 

6 Conclusion 

By combining the classic Bondi model with a more accurate investigation into the asym- 
metric torque, we obtain an analytic model capable of explaining the behaviour of the 
rattleback without numerical calculations. The presentation given here provides a better 
opportunity to understand the dynamics of the rattleback. 
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On the curvature parameters 



The curvature parameters 0, (p and ip originate from Bondi's studies, where their role 
is to simphfy the classification of celt types. This is done by studying the roots of the 
characteristic equation. Bondi defines the parameters as follows: 

e + (j) = p + s 

e^ = ps- g2 (35) 
ip{6 — (j)) = p — s 

According to this definition the parameters are independent of the sign of q or the 
direction of the skew angle. By choosing the sign of q we simultaneously fix the direction 
of the skew angle and the stable direction of spin. 

Choosing q to be positive and solving for the original parameters gives 

p =1(0 + </,) + 1^(0 -</>) 
q=^{e-(P){l-i;^)l (36) 

s = + ^)- 1^1^(6 -cl>). 
By defining p = a/n and Q = g/ av? we can write the characteristic equation as 

(p2 + l)(/)2 + + K + J7A) + f)p + V?v = 0, (37) 
where Bondi defines the new auxiliary parameters as follows 



amx = \[a-m-m-^'')'^ 

am^^ = 1 - l(a + /3 - 27)(^ + <^) + (a - 7)(/3 - l)H - ^(a - /3)(0 - ^ 

1 1 (38) 

= 2^" + ^'^^^ + " ^^'^^ " 2^" " ^'^^^ ~ ^^"^ 
a/30^fi = 2 - (0 + (/.) - (a + /3 - 7)(0 + - 29^) 
a/39(pi^ = (l-6')(l-0) 

The real part of the root cannot change sign without vanishing at some point, which 
enables us to find the critical values separating different celt types. In practice, we end 
up with conditions for the parameters fi and k, which in turn restrict the curvature and 



inertial parameters through equations (38). To draw the diagram mentioned earlier in 
the text, we need to fix some parameters. In the diagram taken from Bondi the values 
of the inertial parameters have been fixed and the division lines given with four different 
values of the skew angle. Celts situated left to the line /i = are unstable in both 
directions, whereas type I celts between the lines = and k = are unidirectional. 
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The more complex type II celts lie right to the line k = 0. As mentioned earlier, the 
diagram shows four different cases with varying skew angle ip. 




The equations used to simplify the analysis are not very interesting as such, but 
they do reveal the conditions of unidirectionality to agree with the ones given by Gar- 
cia&Hubbard. If these requirements are not satisfied, x = the term causing the 
unidirectional behaviour will disappear. The coefficients a, /3, 9, and (p cannot vanish for 
a physical object. 

Although the conditions seem to agree, the curvature parameters involved here have 
been derived in a completely different way. The definition given here stems from their 
convenience in the analysis of roots, whereas in the paper by Garcia&Hubbard these are 
associated with the geometry of the celt. The equations of Bondi are used in the analysis 
together with the ones derived in Garcia&Hubbard, which requires the definitions to be 
equivalent. Establishing this relation is the main topic of this appendix. Although we 
are intrested in confirming the result, the calculation itself is quite illuminating as such. 

To compare the definitions, we need to solve the curvature parameters from the defi- 
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nition given by Bondi: 



= l(p + s) + ^V(p-^)' + V 



(39) 



In the system defined by the axes of inertia, the object has an elhpsoidal form: 



1 - Ipi-f 
2 a 



xy 



2 ^a^ 



This can be expressed in a matrix form as 



12 ll 

2a 2a 

Iq Is 

2a 2a 



The eigenvalues needed in diagonaUzing this matrix are 



a 

Using these, we obtain the eigenvectors 

-9 

1 



a 



and 



1 



(40) 



(41) 



(42) 



(43) 



In the system defined by the eigenvectors, the expression is diagonal. This corresponds 
to a symmetric clHpsoid i.e. the axes we have found arc the axes of curvature. These arc 
orthogonal as required. The other joins the corresponding axis of inertia smoothly as 

— > 0, whereas the other approaches its counterpart when the lower component becomes 
negligible during normalization. 

To find the angle between the axes, we need to normalize the eigenvectors. Let us 
choose the second vector and divide it by its length. Prom the scalar product of the 
normalized eigenvector and its counterpart [0 1]^, we obtain 



cos^ = 



1 + 



(44) 



{l2{p-s)+\^{p-s)'^+Aq^y 

Using the identity cos^ ^ — sin^ ^ = 2 cos^ C ~ 1 results in 



cos^ ^ — sin^ ^ 



(45) 



The definitions of are thus equivalent. 
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For verifying the two other parameters we need the formula for curvature radius 



R 



(46) 



The curvature axes are the eigevectors and according to the general theory, the diagonal 



form contains the eigenvalues. By using formula (46) to the diagonal form, we obtain 
the radii of curvature near the equilibrium to be 



Ri 



a 



and R2 



a 



(47) 



^(p + s) + iV(p-s)2+4g2 Up + '^) - 5\/(p-s)^+ V 

The definition of Garcia&Hubbard thus gives 

= ^ = ^(P+5) + ^\/(p-s)^+ V and(/<= J- = ^{p+s)-^y^{p-sf + Aq^ (48) 

in agreement with the definition by Bondi. 

The stability condition ([s]) requiring the radius of curvature to be larger than the 
equilibrium height of the center of mass can now be understood more properly. This 
requirement is indeed satisfied if 



1 > p, 1 > s and (1 — p){l — s) > 



(49) 



as previously stated. 

Interestingly enough, it seems inevitable that Bondi must have done a similar calcu- 
lation to establish the stability condition. Bondi has probably also solved his definition 
for the curvature parameters. From these two results the connection of the parameters 
to the celt shape is quite easy to see, but Bondi does not mention anything about their 
geometrical meaning. This is in retrospect easy to see, but fairly difficult to see directly 
from the equations without precognition, which may have caused it to remain unnoticed. 
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